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Abstract 

It is shown how to model weakly dissipative free-surface flows using the classical potential flow approach. The 
Helmholtz-Leray decomposition is applied to the linearized 3D Navier-Stokes equations. The governing equations 
are treated using Fourier-Laplace transforms. We show how to express the vortical component of the velocity only 
in terms of the potential and free-surface elevation. A new predominant nonlocal viscous term is derived in the 
bottom kinematic boundary condition. The resulting formulation is simple and does not involve any correction 
procedure as in previous viscous potential flow theories [7]. Corresponding long wave model equations are derived. 

Resume 

Ecoulements viscopotentiels avec surface libre en profondeur finie. Nous montrons comment modeliser 
les ecoulements a surface libre faiblement dissipatifs en utilisant une approche potentielle. La decomposition de 
Helmholtz-Leray est appliquee aux equations de Navier-Stokes linearisees. Le probleme est etudie au moyen de la 
transformee de Fourier-Laplace. Nous montrons comment exprimer la partie rotationnelle de la vitesse en fonction 
du potentiel des vitesses et de I'elevation de la surface libre. Un nouveau terme nonlocal preponderant apparait 
dans la condition cinematique au fond. La formulation finale est simple et ne requiert pas de corrections ulterieures 
comme dans [7]. Un modele d'ondes longues est obtenu a partir de ces nouvelles equations. 

Key words: potential flow, free-surface flow, viscosity, dissipation, water waves, wave damping 

Mots-cles : ecoulement potentiel, ecoulement avec surface libre, viscosite, dissipation, vagues, attenuation de vague 



Version frangaise abregee 



Recemment, Dias et al. [3] ont donne une nouvelle formulation du probleme des vagues en profondeur 
infinie en presence de faible dissipation. L'idee principale est de partir des equations de Navier-Stokes 
linearisees (1) oii v represente la viscosite cinematique, v le vecteur vitesse, g I'acceleration due a la 
gravite, p la pression, p la densite, t le temps, et d'utiliser la decomposition de Helmholtz-Leray (2) pour 
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la Vitesse. Ensuite I'on exprime la composante rotationnelle 'ip de la vitesse en fonction de la composante 

potentielle (j> et de la deformation rj de la surface libre. Nous presentons ici une extension au cas d'une 
couche de fluide de profondeur finie. Dans ce cas-la il faut tenir compte des frottements au fond. 

En utilisant la transformee de Fourier-Laplace, on ecrit la solution generale pour (p et tp. Ensuite Ton 
utilise la condition dynamique sur la surface libre. La nuUite des deux composantes tangentielles de la 
contrainte donne les equations (3) et (4), qui peuvent etre recombinees pour donner I'equation (5). Cette 
meme equation, combinee avec la condition cinematique sur la surface libre (6), donne un lien entre 
r] et xjj (7) et un autre lien entre 77 et (8). Le remplacemcnt dc la composante rotationnelle dans la 
condition cinematique a la surface libre donne la premiere equation modifiee (9). De memo I'equilibre de 
la contrainte normalc au niveau de la surface libre donne I'equation de Bernoulli modifiee (10). 

La scconde partie de la note consistc a incorporcr une correction dc type couclic limitc dans la condition 
cinematique au fond. L'idce est d'introduirc Tcpaisscur dc la couche limitc proportionncUe h. 5 — y^, la 
coordonnee dc couche limite C, = {z + h)/S, puis d'effectuer un developpement du potcntiel et de la vitesse 
rotationnelle en fonction du petit paramctrc S. On obtient ensuite une hierarchic de problemes a resoudre. 
Ce qui nous interesse, c'est la correction a apporter a la condition cinematique au fond (13). Aprcs quelques 
calculs, on obtient la correction presentee dans I'equation (14). Finalcment nous proposons Ic nouveau 
systeme d'equations (15)-(18) pour I'etude des vagues avec faible viscosite en profondeur finie. Notons 
que les termes nonlineaires ont ete ajoutes de fagon heuristique. A partir de ce systeme, nous pouvons 
deriver le systeme de Boussinesq faiblement dissipatif donne par les equations (19)-(20). 



1. Introduction 



The effects of viscosity on gravity waves have been addressed since the end of the nineteenth century 
in the context of the linearized Navier-Stokes (NS) equations. It is well-known that Lamb [9] studied this 
question in the case of oscillatory waves on deep water. What is less known is that Boussinesq studied 
this effect as well [1]. In this particular case they both showed that 



ot 



where a denotes the wave amplitude, v the kinematic viscosity of the fluid and k the wavenumber of the 
decaying wave. This equation leads to the classical law for viscous decay, namely a{t) = aoe"^"^ *. 

In this work we keep the features of undamped free-surface flows while adding dissipative effects. The 
classical theory of viscous potential flows [7] is based on pressure and boundary conditions corrections 
due to the presence of viscous stresses. We present here a novel approach. 

Currently, potential flows with ad-hoc dissipative terms are used for example in direct numeric;al sim- 
ulations of weak turbulence of gravity waves [5,6,11]. There have also been several attempts to introduce 
dissipative effects into long wave modelling [10,4,2]. 

The present article is a direct continuation of the recent study [3] . In that work the authors considered 
two-dimensional (2D) periodic waves in inflnite depth, while in the present study we remove these two 
hypotheses and all the computations are done in 3D. This point is important since the vorticity structure 
is more complicated in 3D. In other words we consider a general wavetrain on the free surface of a fluid 
layer of flnite depth. As a result we obtain a qualitatively different formulation which contains a nonlocal 
term in the bottom kinematic condition. The inclusion of this term is natural since it represents the 
correction to potential flow due to the presence of a boundary layer. Moreover, this term is predominant 
since its magnitude scales with 0{^), while other terms in the free-surface boundary conditions are of 
order 0{v). Other researchers have obtained nonlocal corrections but they differ from ours [8]. 
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2. Derivation 



Consider the hnearized 3D incompressible NS equations describing free-surface flows in a fluid layer of 
uniform depth h: 

^ = — Vp + uAv+g, V-v = 0, (1) 

with V the velocity vector, p the pressure, p the fluid density and g the acceleration due to gravity. We 
represent v = [u, v, w) in the form of the Helmholtz-Leray decomposition: 

v = V^ + V xtp, V = (V'l,V'2,^3). (2) 

After substitution of the decomposition (2) into (1), one notices that the equations are verified provided 
that the functions (j) and ip satisfy the following equations: 

p at 

Next we discuss the boundary conditions. We assume that the velocity field satisfies the conventional 
no-slip condition at the bottom v\^^_f^ = 0, while at the free surface we have the usual kinematic 
condition rit = w and three dynamic conditions [ct • n] =0, where a is the stress tensor, [/] denotes the 
jump of a function / across the free surface, and the normal vector n equals (0, 0, 1) due to linearization. 

Using Fourier-Laplace transforms, which we denote by £^ = CoJ^, f{x, t) /(fe, s), k = (fc^;, ky) we 
can determine the structure of the unknown functions (/), xp in the transform space. We assume that all 
the functions involved in the present computation satisfy the necessary regularity requirements and have 
sufficient decay at infinity so that the integral transforms can be applied. The solution for ^ is obtained 
from the transformed continuity equation A(p ~ — ^ ^zz — \kf 4^ = and ip from the corresponding 
transformed equation ip^ = uAip sip = v(ipzz ~ 1^1^ V')- 

4> = <^+(fe,s)el'=l^ + (fe.s)e"l'°l", i^i = V'io(fe,s)(el™l^ + Q(fe, 5)6-1"!^), 

where := \kf + s/u and (p^, (p^, ip^, C := (Ci,C2,C3) are unknown functions of the transform 
parameters (fe, s), determined by the initial and appropriate boundary conditions. 

There are three dynamic conditions on the free surface. Let us use first those related to the tangential 
stresses (the third one will be used later), where p, = pv: 

( dw du\ f dw dv\ 

Substituting decomposition (2) into these two identities yields 

^ 92^-2 d^iPi 92^3 d^^2 



dxdz dx^ dxdy dydz dz"^ 



0, = 0, (3) 



92^2 Q2^^ Q2^^ Q2^, 



3=0, = 0. (4) 



dydz dxdy dy^ dz'^ dxdz 
The next step consists in taking the Fourier-Laplace transform to these relations. We do not give here 
the explicit expressions since this operation is straightforward. The combination {—ikx){S) + {—iky){4:) 
gives the important relation 

ikyi^wil + Ci) - ik j2o{l + C2) = p^^i'Po - V>o)- (5) 

m'^ -I- fc 
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Let us turn to the free-surface kinematic condition ^ = w = ^ + — ^ = 0. In transform 
space it becomes 

si) = \k\ {ip^ - iPq) + ikyipio{l + Ci) - ifcxf/'2o(l + C'2). (6) 
Equations (5) and (6) can be rewritten as 

mi^^fj, (7) 

ikyi)w{l + Ci) - ifc^V'2o(l + C2) = -1v r]. (8) 

Using (8) one can replace the rotational part in the kinematic free-surface condition: 

% =0. +>C^'[-2i^|fc|'r)] =</'^ + 2i^A77. (9) 

In order to account for the presence of viscous stresses, we have to modify the dynamic free-surface 
condition as well. This is done using the balance of normal stresses at the free surface: 

dw fd'^cf) d'^ijj2 \ 

= at z = =^ p - Po = ^pv— = 2pj^l ^-^ + -r-r- - ^-r- 1. 

az \oz'^ axaz avaz/ 



Using (8) one can show that gjgt — dv^z — ^iy^ \ Bernoulli's equation becomes 



dz \dz'^ dxdz dydz' 

dxdz 

<i)t+g'n + 2v(l),, + 0{vi) = Q. (10) 

Since we only consider weak dissipation [v ^ 10^^ — lO^"* m^/s). wc neglect terms of order o{u). 

The second step in our derivation consists in introducing a boundary layer correction at the bottom. 
Obviously, this was not done in the previous study [3], since the derivation dealt with the infinite depth 
case;. In order to include this modification, we consider a semi-infinite fluid layer as it is usually done in 
boundary layer theory. The fluid occupies the domain z > —h. In this derivation we use the pure Leray 
decomposition of the velocity field v = V(p + u together with the divergence-free constraint V • ti = 0. 
Expecting that the rotational part u varies rapidly in a distance S = y/V, ^ we introduce the boundary- 
layer coordinate ( = {z + h)/S, so that u = u{x,() = {ux,Uz){x,(). The solid boundary is given by 
C = 0, and the potential part of the flow is not subject to this change of variables. With the new scaling, 
the divergence-free condition becomes 

^ + 6V^-u^ = 0. (11) 

As done in [10], we expand the imknown fmictions in powers of the small parameter 6: 
(j) = (j)o{x, z, t) + 5(j)i{x,z,t) + u = Qaix, (, t) + 5qj^{x, C,t) + .... 

Substituting the expansion for u into (11) gives the following relations: 

where Qq denotes the first two components of the vector (/q corresponding to the horizontal coordinates 
x. Recall that we would like to determine the correction to the bottom boundary condition 

(t>z = - Uzl^^g = - (go, +'5'7iJ|^=o +o((5). (13) 

So we only need to compute qo^ and qi^ at the bottom C = 0- 



^ Of course we should nondimensionalize all quantities in order to define small numbers. One would find that S is in fact 
equivalent to VRe~^ x L, where Re is the Reynolds number and L a typical length. 



Using the same asymptotic considerations as above, we can write down the following sequence of 
problems: 



A0o = O, ^ 



oz 



1. 



together with the radiation condition q — > as C ^ co- 

This sequence of linear problems can be solved using Fourier transforms. In Fourier space one finds 
immediately that (^o(i, z, k) = (po{t, k) (el*"!^ + e"!**!^) . Since we know ^q, we can determine the rotational 
component q^. 

Analytical solutions to the equation dq^^/dt = d'^qQ^/d('^ are well-known. If we assume that initially 
the flow is potential and the boundary condition is q^^ = ik^o{z = —h;k), the solution is 

t 



Iff - 
—= / -3-e ^i^--)ik(t>o{T,z = -h,k)dT. 





Let us now integrate the second equation in (12) from to oo, using the appropriate decay at infinity: 

oo oo t 

9iJc=o = -/ ^''•^o. '^C=^ J J -—^e-^\k\''MT,z = -h,k)dTdC. 



One can interchange integral signs and evaluate the inner integral on ^ to obtain: qi^ = J ^''^ 4>o(T^z—-h,k) 
Hence, the bottom boundary condition becomes, at order 6, 



z=—h 



dT=J- / - dT = -W- / dr. (14) 



\Jt — T V I" 7 \/t — T V J \/t — T 



One recognizes on the right-hand side a half-order integral operator. Summarizing the developments made 
above and generalizing our equations by including nonlinear terms (this is a conjecture at this stage), we 
obtain a new set of viscous potential free-surface flow equations: 

A(f> = 0, {x,z) en = R'^ x[-h,ri\ (15) 

m + Vr/- V0 = 0^ -|-2!/A?7, z = r] (16) 

(pt + ^\'7(pf+9V = -'^^(t>zz, z = V (17) 



f 4>zz 



dr, z = -h. (18) 



Using this weakly damped potential flow formulation one can derive the following system of Boussinesq 
equations with horizontal velocity Uh defined at the depth z$ = —Oh, < ^ < 1: 



t 

Vt + V-iih + v)uh) + h^(^^-e+^^V\V-Uh) = 2uATj+^J^=^dT, (19) 



Uht + ^V|uh|2 + gVv - h?e (1 - ^ ) V(V • Uht) = IvAuh- (20) 
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3. Conclusion 



In the present paper we have shown how to express the rotational component of the velocity field in 
terms of the potential part of Helmholtz-Leray decomposition. This expression contains differential and 
integral operators. Obviously, this analysis is only linear. In future work wc will try to extend the present 
derivation to the nonlinear case. A long wave approximation was derived from this new potential flow 
formulation. 

It is interesting to note that dissipativc terms of this form have been used to verify the theory of weak 
turbulence of surface gravity waves in deep water [5,6,11]. They were added without justification to model 
dissipation at small scales. Note that a good qualitative agreement was obtained between the Kolmogorov 
spectrum predicted by weak turbulence theory and the results of DNS. Hence, the present work can be 
considered as an attempt to justify the inclusion of these terms. 

Our final remark concerns the nonlocal term in the kinematic bottom boundary condition. This term 
can be also considered as a boundary layer correction at the bottom. In modelling viscous effects this 
term plays the main role, since its magnitude is C'(v^). Of course, the numerical implementation of this 
term is another matter. 

What is the value of v to be taken in numerical simulations? There is surprisingly little published 
information of this subject. What is clear is that the molecular diffusion is too small to model true 
viscous damping and one should rather consider the eddy viscosity parameter. 
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